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ABSTRACT 

We present a theoretical study of radiative transfer in strongly magnetized electron-ion plasmas, 
focusing on the effect of vacuum polarization due to quantum electrodynamics. This study is directly 
relevant to thermal radiation from the surfaces of highly magnetized neutron stars, which have been 
detected in recent years. Strong-held vacuum polarization modihes the photon propagation modes in 
the plasma, and induces a “vacuum resonance” at which a polarized X-ray photon propagating outward 
in the neutron star atmosphere can convert from a low-opacity mode to a high-opacity mode and vice 
versa. The effectiveness of this mode conversion depends on the photon energy and the atmosphere 
density gradient. For a wide range of held strengths, 7 x 10^^ ^ B ^ a few x 10^® G, the vacuum 
resonance lies between the photospheres of the two photon modes, and the emergent radiation spectrum 
from the neutron star is signihcantly modihed by the vacuum resonance. (For lower held strengths, only 
the polarization spectrum is affected.) Under certain conditions, which depend on the held strength, 
photon energy and propagation direction, the vacuum resonance is accompanied by the phenomenon 
of mode collapse (at which the two photon modes become degenerate) and the breakdown of Faraday 
depolarization. Thus, the widely used description of radiative transfer based on photon modes is not 
adequate to treat the vacuum polarization effect rigorously. We study the evolution of polarized X-rays 
across the vacuum resonance and derive the transfer equation for the photon intensity matrix (Stokes 
parameters), taking into account the effect of birefringence of the plasma-vacuum medium, free-free 
absorption, and scatterings by electrons and ions. 

Subject headings: magnetic helds - radiative transfer - stars: neutron - stars: atmospheres - X-rays: 
stars 


1. INTRODUCTION 

It has long been known that in a strong magnetic 
held even the vacuum has a nontrivial dielectric prop¬ 
erty; this is the result of vacuum polarization, a genuine 
quantum electrodynamics effect (Heisenberg & Euler 1936; 
Schwinger 1951; Adler 1971; Tsai & Erber 1975; Heyl 
& Hernquist 1997). This ehect is usually negligible for 
B ^ Bq, but becomes increasingly important when B ex¬ 
ceeds Bq, where Bq = m‘lc^/{eh) = 4.414 x 10^^ G is 
the magnetic held at which the electron cyclotron energy 
hujBe = heB/ (mec) equals meC^. A number of studies 
in the late 1970s showed that vacuum polarization can af¬ 
fect radiative transfer in magnetized plasmas by modifying 
photon polarization modes and inducing features in the 
radiative opacities (Gnedin, Pavlov & Shibanov 1978a,b; 
Meszaros & Ventura 1978,1979; Pavlov & Shibanov 1979; 
Ventura, Nagel & Meszaros 1979; see Pavlov & Gnedin 
1984 and Meszaros 1992 for review). These earlier studies 
were mostly applied to radio pulsars and accreting X-ray 
binaries, for which B ^ Bq. Not surprisingly, the ef¬ 
fect of vacuum polarization in these systems appears to be 
relatively minor and, as far as we know, no concrete ob¬ 
servational manifestation of vacuum polarization has been 
identihed. 

In the last few years, soft gamma-ray repeaters (SGRs) 
and anomalous X-ray pulsars (AXPs) have been identihed 
as a potentially new class of neutron stars (NSs) endowed 
with superstrong {B ^ 10^"^ G) magnetic helds, the so- 
called “magnetars” (see Thompson & Duncan 1996; Hur¬ 
ley 2000; Thompson 2001; Israel et al 2001). This has 


stimulated new studies on radiative processes in the su¬ 
perstrong magnetic held regime. Of particular relevance 
to the present paper is the detection of thermal X-rays 
from the surface of isolated magnetic NSs (see Pavlov et 
al. 2002 for a recent review), including hve magnetar can¬ 
didates (4 AXPs and one SGR; see, e.g., Patel et al. 2001; 
Juett et al. 2002; Tiengo et al. 2002). The thermal emis¬ 
sion can potentially provide useful constraints on the in¬ 
terior physics, surface magnetic held and composition of 
NSs. Several recent works have begun to model the at¬ 
mosphere emission from magnetars, including the ehect 
of vacuum polarization (Bezchastnov et al. 1996; Bulik & 
Miller 1997; Ozel 2001; Ho & Lai 2001, 2003; Lai & Ho 
2002a; a review/comment of these works can be found in 
Ho & Lai 2003). 

In a highly magnetized plasma which characterizes NS 
atmospheres, there are two photon polarization modes, 
and they have very diherent radiative opacities. Previ¬ 
ous theoretical modeling of magnetic NS atmospheres (see 
Ho & Lai 2003 and references therein; see also Pavlov et 
al. 1995; Zavlin & Pavlov 2002 for reviews) are based on 
solving the transfer equations for the two photon modes. 
As we discuss and quantify in §2 and §3 below, this modal 
description of the radiative transport is, in general, inad¬ 
equate to handle the vacuum polarization ehect near the 
“vacuum resonance” where the modes can mix and even 
“collapse”. The main purpose of this paper is to clarify 
these issues and to derive general transfer equations that 
can be used to study the transport of polarized radiation 
in the superstrong held regime. 
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Section 2 contains a physical discussion of the vacuum 
polarization effect on radiative transfer in magnetized NS 
atmospheres. In §3, we derive the basic properties of pho¬ 
ton modes in magnetized plasmas and quantify the limi¬ 
tation of the modal description of radiative transport. We 
then study in §4 and §5 the propagation of polarized pho¬ 
tons in an inhomogeneous plasma. We derive in §6 the 
transfer equations for the photon intensity matrix and con¬ 
clude in §7. 

2. VACUUM POLARIZATION EFFECT ON RADIATIVE 
TRANSFER AND ATMOSPHERE SPECTRA OF 
MAGNETARS: A PHYSICAL DISCUSSION 

To motivate our quantitative study in §3-§6, we present 
a physical discussion of the effect of vacuum polarization 
on radiative transfer in strong magnetic fields and on the 
thermal spectra of magnetars (see Lai & Ho 2002a; Ho & 
Lai 2003; hereafter LH02 and HL03). This also serves to 
clarify some of the confusions that have appeared recently 
(see Lai & Ho 2002b). 

In a highly magnetized NS atmosphere, both the plasma 
and vacuum polarization contribute to the dielectric prop¬ 
erty of the medium. A “vacuum resonance” arises when 
these two contributions “compensate” each other (Gnedin 
et al. 1978a,b; Pavlov & Shibanov 1979; Meszaros & Ven¬ 
tura 1979; Ventura et al. 1979). For a photon of energy 
E, the vacuum resonance occurs at the density 

pv^QMAY-^Bl^Elr^ gcur\ ( 1 ) 

where Te is the electron fraction, Ei = E/{1 keV), B = 
10^‘*i?i4 G is the magnetic field strength, and / = f{B) 
is a slowly varying function of B and is of order unity 
(LH02 and HL03). For p > pv (where the plasma ef¬ 
fect dominates the dielectric tensor) and p < pv (where 
vacuum polarization dominates), the photon modes (for 
E much smaller than the electron cyclotron energy Ebb) 
are almost linearly polarized (see Fig. 3 below): the ex¬ 
traordinary mode (X-mode) has its electric field vector E 
perpendicular to the k-B plane, while the ordinary mode 
(0-mode) is polarized along the k-B plane (where k spec¬ 
ifies the direction of photon propagation, B is the unit 
vector along the magnetic field). Near p = pv, however, 
the normal modes become circularly polarized as a result 
of the “cancellation” of the plasma and vacuum effects — 
both effects tend to make the mode linearly polarized, but 
in mutually orthogonal directions. When a photon propa¬ 
gates in an inhomogeneous medium, its polarization state 
will evolve adiabatically if the density variation is suffi¬ 
ciently gentle. Thus, a X-mode (0-mode) photon will be 
converted into a 0-mode (X-mode) as it traverses the vac¬ 
uum resonance (see Fig. 1). This resonant mode conver¬ 
sion is analogous to the MSW effect of neutrino oscillation 
(e.g., Haxton 1995). For this conversion to be effective, 
the adiabatic condition must be satisfied: 

E^E,d{B,eB,Hp) 


where 9b is the angle between k and B, m = {EBi/E)‘^, 
Ebi is the ion cyclotron energy, and EIp = \dz/d\n.p\ 
is the density scale height (evaluated a,t p = pv) along 
the ray. For an ionized Hydrogen atmosphere, EIp ~ 
2kT/{mpg cos 9) = 1.65 Te /((714 cos 0) cm, where T = 
10® Tq K is the temperature, g = 10 ^^( 7 i 4 cm s“^ is the 
gravitational acceleration, and 9 is the angle between the 
ray and the surface normal. We note that in this paper [as 
in our previous papers (LH02, HL03), and in several other 
references such as Meszaros & Ventura (1979), Meszaros 
(1992, p. 94), etc.], we refer to the 0-mode as the mode 
with \K\ = \Ex/Ey\ ^ 1 and the X-mode as the mode 
with jATj <C 1 (here the z-axis is along k and the y-axis is 
in the direction of B x k), thus the name “mode conver¬ 
sion” . Alternatively, one can define modes with definite 
helicity (the sign of AT = —iE^/Ey) such that K changes 
continuously as p changes; we call these plus-mode and 
minus-mode (see §3). Thus we may also say that in the 
adiabatic limit, the photon will remain in the same plus or 
minus branch, but the character of the mode is changed 
across the vacuum resonance. Indeed, in the literature 
on radio wave propagation in plasmas (e.g., Budden 1961; 
Zheleznyakov et al. 1983), the nonadiabatic case, in which 
the photon state jumps across the continuous curves, is 
referred to as “linear mode coupling”. It is important to 
note that the “mode conversion” effect discussed here is 
not a matter of semantics. The key point is that in the 
adiabatic limit, the photon polarization ellipse changes its 
orientation across the vacuum resonance, and therefore the 
photon opacity changes significantly We prefer to use 
the term “mode conversion” since it captures this essential 
feature important for radiative transfer, and it is analogous 
to the MSW effect and other quantum mechanics problems 
involving adiabatic evolution of quantum states (Landau 
1932; Zener 1932). 

Because the two photon modes have vastly different 
opacities, the vacuum polarization-induced mode conver¬ 
sion can significantly affect radiative transfer in mag- 
netar atmospheres. The main effect of vacuum po¬ 
larization on the atmosphere spectrum can be under¬ 
stood as follows. When the vacuum polarization ef¬ 
fect is neglected, the decoupling densities of the 0-mode 
and X-mode photons (i.e., the densities of their respec¬ 
tive photospheres) are approximately given by (for Hy¬ 
drogen plasma) po — g cm“® and 

px — 486Tg~^^^i?y^Hi4G“^/^ g cm“® (see LH02), where 
G = 1 —The vacuum resonance lies between these 
two photospheres when po < Pv < Px, i-e. 

Bi<B < Bh, (3) 

where 

Bi ~ 6.6 X 10^® /G, (4) 

Bh ~ 5.1 X 10^® f T~^I^E~'^^'^G-^I'^ G. (5) 

When this condition is satisfied, the effective decoupling 
depths of the photons are changed^. Indeed, we see from 


= 1.49 (/ tan6»B|l - 

^ The O-mode has a significant component of its E field along B (for most directions of propagation except when k is nearly parallel to B), 
and therefore the O-mode opacity is close to the B = 0 value, while the X-mode opacity is much smaller. 

^ For B > Bh the vacuum resonance lies deeper than the photospheres of both modes, while for B < Bi the resonance lies outside both 
photospheres. In both cases, the effect of vacuum polarization on the radiative spectrum is expected to be small, although in the latter case 
(B < Bi), the polarization of the emitted photons will be modified by the vacuum resonance. 
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Fig. 2 that mode conversion makes the effective decoupling 
density of X-mode photons (which carry the bulk of the 
thermal energy) smaller, thereby depleting the high-energy 
tail of the spectrum and making the spectrum closer to 
black-body (although the spectrum is still harder than 
black-body because of nongrey opacities) This expec¬ 
tation is borne out in self-consistent atmosphere modeling 
presented in HL03. Another important effect of vacuum 
polarization on the spectrum, first noted in HL03, is the 
suppression of proton cyclotron lines (and maybe other 
spectral lines; see also Ho et al. 2003). The physical ori¬ 
gin for such line suppression is related to the depletion of 
continuum flux, which makes the decoupling depths inside 
and outside the line similar. HL03 suggests that the ab¬ 
sence of lines in the observed spectra of several AXPs may 
be an indication of the vacuum polarization effect at work 
in these systems. 

Our previous study (LH02) on the vacuum-induced 
mode conversion did not explicitly take into account of 
the effect of dissipation on the mode structure. Although 
this dissipative effect is small under many situations, near 
the vacuum resonance and for some directions of propaga¬ 
tion, the two photon modes can “collapse”, i.e., they be¬ 
come identical and hence nonorthogonal (Soffel et al. 1983; 
see also Pavlov & Shibanov 1979). The analysis of LH02 
breaks down near these “mode collapse” points (see §3). 
The question therefore arises as to how the polarization 
state evolves as the photon travels near the mode collapse 
point. Also, as mentioned in §1, previous studies of NS at¬ 
mospheres with strong magnetic fields rely on the modal 
description of the radiative transport. This is valid only in 
the limit of large Faraday depolarization (Gnedin & Pavlov 
1974), which is not always satisfied near the vacuum res¬ 
onance especially for superstrong field strengths (see §3). 
More importantly, the transfer equations based on nor¬ 
mal modes cannot handle the cases in which partial mode 
coupling (conversion) occurs across the vacuum resonance 
(i.e., when the adiabatic condition is neither strongly sat¬ 
isfied or violated). Clearly, to properly account for the 
effects of mode collapse, breakdown of Faraday depolar¬ 
ization, and mode conversion associated with vacuum po¬ 
larization, one must go beyond the modal description of 
the radiation field by formulating and solving the trans¬ 
fer equation in terms of the photon intensity matrix (or 
Stokes parameters) and including the birefringence of the 
plasma-vacuum medium (see §5 and § 6 ). 

3. PHOTON POLARIZATION MODES, VACUUM 
RESONANCE, MODE COLLAPSE, AND FARADAY 
DEPOLARIZATION 

The property of photon modes in a magnetized electron- 
ion plasma including vacuum polarization was studied in 
LH02 and HL03 (see also Meszaros 1992 and references 
therein for earlier works in which electron plasmas with 
B <C Bq were studied). In these papers, however, the 
damping terms in the dielectric tensor were largely ne¬ 


glected (but see §2.5 of HL03). Including these damp¬ 
ing terms can (for certain parameter regimes) modify the 
mode property near the vacuum resonance and give rise 
to the phenomenon of mode collapse. 

Following Ginzburg (1970), we consider a cold, magne¬ 
tized plasma composed of electrons and ions (with charge, 
mass and number density given by —e, me, rze and 
Ze, M = Amp, rii = nejZ, respectively; here Z is the 
charge number and A is the mass number of the ion). 
The electrons and ions are coupled by collisions (with the 
collision frequency Pei)- We generalize Ginzburg’s result 
by including the radiative dampings of electrons and ions, 
with the damping frequencies Vre and v^i, respectively. In 
the coordinate system XYZ with B along Z, the plasma 
contribution to the dielectric tensor is given by ^ 



where 

Veil + i"fri) + Vi{l + i^re) 

e±g~l -(7) 

(1 -I- i^re ± uj )(1 -I- I'^ri T u/ ) -|- 


ry ~ 1 — 


1 + z(7ei + 7re) 1 + z(7ei + 7ri) ' 


In eqs. (7) and ( 8 ), we have defined the dimensionless 
quantities 

Ue = ^, Ve=‘^, = (9) 

where uiBe = eB/iniec) is the electron cyclotron fre¬ 
quency, iOBi = ZeB/iMc) = iZme/Amp)ujBe is the ion 
cyclotron frequency, tOpe = {YnUe^jviieY^'^ is the elec¬ 
tron plasma frequency, and Wp, = (dTrrziZ^e^/M)^/^ = 
{ZnielAmpY/'^ijjpe is the ion plasma frequency. The di¬ 
mensionless damping rates 7 ei = Ire = VrejoJ, and 

7 ri = Vrijuj are given by 





z'^me'^ ( 

lei 

~ \ 


= 9.2 X 10“ 



Ire 

3m eC^ 


Z'^me 

Iri 

Amp 


A T^'^El 


where is the Gaunt factor. In eq. (7), 

lfi=lei ■2’"^)+me/(Amp) . (13) 

For Hydrogen plasmas, 7 ^ ~ 7 ei. 

Including vacuum polarization, the dielectric tensor e 
and inverse permeability tensor = fi can be written 


^ Even when mode conversion is neglected, the X-mode decoupling depth can still be affected by vacuum polarization. This is because the 
X-mode opacity exhibits a spike feature near the resonance, and the optical depth across the resonance region can be significant; see LH02. 

^ Unlike the expression given in Ginzburg (1970), equations (7) and (8) are accurate as long as 7 ei, 7re, 7ri ^ 1; these accurate expressions 
are necessary in order to recover the correct scattering cross-sections. Also note that when the damping terms are neglected, eqs. (7)-(8) agree 
with eqs. (2.2)-(2.4) of HL03. When the damping terms are included, the latter equations are incorrect. The reason is that in the presence of 
electron-ion collisions, the electron and ion cannot be treated as independent particles. The difference between these two sets of equations is 
appreciable only for u) ^ u)Bi- See Potekhin Chabrier (2002) for a discussion. 
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as e = + Ae^'"'>, fj, = 1 + (where I is the unit 

tensor). The vacuum corrections are 

= (a - 1)1 + gBB, (14) 

A/i(’') = (a - 1)1 + mBB, (15) 

where a, q, and m are functions of B (the explicit expres¬ 
sions are given in LH02 and HL03). Including the contri¬ 
bution from vacuum polarization, the dielectric tensor can 
still be written in the form of eq. (6), except 

£^£' = £- 1 - 0 — 1 , 'q^ri'=ri + a-\-q—1. ( 16 ) 


With the dielectric tensor and permeability tensor 
known, the photon modes can be determined in a straight¬ 
forward manner. In the xyz coordinates with k along the 
2 -axis and lies in the X-Z plane (such that B x k = sin 6*^ y, 
where 9b is the angle between k and B), we write the elec¬ 
tric field of the mode as E± oc {iK±, l,iKz±)^ where the 
ellipticity K± = —iE^jEy is given by 


= /3 ± V /32 + r, (17) 


with r = 1 -|- (raja) sin^ 9b — 1, and the (complex) polar¬ 
ization parameter (3 is 


/?=- 


£'^ — — e'r]' (1 -I- m/a) sin^ 9b 


2gr]' 


cos 9 b 


The refractive index n± is given by 


2 9V 
n± = 


1 


0633 \g K± 


cos 9 b 


(18) 


(19) 


where £33 = e' sin^ 9b +g' cos^ 9b- 

The polarization parameter {3 directly determines the 
characteristics of photon normal modes in the medium. 
For Ve 1, the real part of (3 can be written as Re(/3) = 
/3oPv: where 


and 


Po 


1/2 ■ 2n 
Ue sm 9 b 

2 COS 9b 


(1 - Ui) , 


Pv — 1 -|- 


{q -I- m) (1 - Ue) 

UeVe 


( 20 ) 

( 21 ) 


The condition Re(/3) = 0 selects out three critical photon 
energies Ebb = ftzoBe, Ebi = hiosi and Ey (for a given 
density), where the latter is the vacuum resonance energy 


EyZ, 


HiOpe 

y/q + m 


( 22 ) 


We shall concentrate on this vacuum resonance in the re¬ 
mainder of this paper. Since Ey depends on density, a 
photon with a given energy E traveling in an inhomo¬ 
geneous medium encounters the vacuum resonance at the 
density py, given by eq. (1). Figure 3 show some examples 
of the mode properties near the vacuum resonance. 

Including the damping terms in the dielectric tensor 
gives rise to the phenomenon of mode collapse (see Soffel 
et al. 1983 and references therein). This occurs when the 
two polarization modes become identical, i.e., Ky = K- or 
n_|_ = n_, which requires P = ±i^/r « ±z. Using eq. (18), 
we find, for t6e ^ 1, We ^ 1, and neglecting ions {ui <C 1), 


a 


1/2 ■ 2 
Ue sm 


9b 


{Pv - he) 


where Pv zz 1 — (q + m)/ve = 1 — {E/Ey)"^ and 7 e = 
lei + Ire- So the Condition P = ±i selects out a crit¬ 
ical angle 9b = ^coii at which mode collapse occurs: 

sin^ 0coii/cos6lcoii = 2/{ul^^"fe)- When the ion effect is 
included (i.e., when Ui ^ 1), no simple expression for the 
complex P can be obtained, and 0coii must be calculated 
numerically (see Fig. 4). Figure 4 shows that for “ordi¬ 
nary” field strengths {Bn ^ 0.1), dcoii is very close to 
90° for most photon energies (in agreement with Soffel et 
al. 1983); for superstrong field strengths {Bn ^ 1), how¬ 
ever, 0coii can be much smaller than 90°. Figure 4 also 
shows that at i? = Ebi (the ion cyclotron energy), the 
mode collapse point is always very close io 9 b = 90°. 

More generally, the modal description of radiative trans¬ 
fer is valid only in the limit of large Faraday depolariza¬ 
tion, when the phase shift between the two modes over 
a mean free path is significant, i.e., when the condition 
|Re(n+ — n_)| » |Im(n+ -I- n_)| is satisfied (Gnedin & 
Pavlov 1974). Consider first the situation where the ion 
effect is neglected (i.e., Ui <C 1). For Ue ^ 1, we find 

Re(n+ - nh) « cosdsRe , (24) 

tie 


Im(n+ -I- n?.) « 7 eUe 


Sin 


(1 


COS' 


' 9b) 


(25) 


where P is given by eq. (23). For |/3v| ^ a few x 7 e (i.e., 
for E slightly away from Ey), we have |Re(n^ — rP_)\ ~ 

2veUe ' max{cos 9B, Ue Py sin 9 b) , thus the condition 
for Faraday depolarization is satisfied for all 9b- On the 
other hand, for |/3v| 'C 7e (i-C., for E extremely close to 
Ey), the condition |Re(n+ — n_)| > |Im(n+ -|- n_)| re¬ 
quires ( 7 euy^/ 2 )| sin^ 9b/cos9b\ = |Im(/3)| < 1; this con¬ 
dition translates to 9b P: 6*coii or 0 b > tt — 0coii- In other 
words, the breakdown of Faraday depolarization at the 
vacuum resonance (|/3v| 7e) is restricted to an angular 

region of 0b around 7r/2. For a given B and 0b, to satisfy 
Faraday depolarization requires E :§> Eye,^{B,9B), where 
EFar{B,9B) is obtained by setting Im(/3) = 1. This result 
can be generalized to the case where the ion effect is not 
negligible (i.e., Ui P 1), although the simple expressions 
(24) -(25) are no longer valid. In general, while mode col¬ 
lapse at the vacuum resonance(/3v' = 0) occurs only at two 
angles 0b = 6*coii and tt — 0coii, the breakdown of Faraday 
depolarization occurs for 0con 9b P, tt — 0coii- To put it 
another way, for a given B and 0b, there is a range of ener¬ 
gies for which Faraday depolarization breaks down; we see 
from Fig. 5 that this range becomes larger as B increases. 

The above analysis shows that in the superstrong field 
regime {Bn ^1), mode collapse and breakdown of Fara¬ 
day depolarization at vacuum resonance can occur for a 
wide range of photon energies and propagation directions. 
Therefore, a rigorous treatment of radiative transfer near 
the vacuum resonance, in general, requires solving the 
transport equations for the four photon intensity matrix, 
or Stokes parameters. 


4. EQUATIONS FOR THE EVOLUTION OF 
ELECTROMAGNETIC WAVES AND MODE AMPLITUDES 

An electromagnetic (EM) wave with a given frequency 
(jj satisfies the wave equation 

^2 

V X (^ • V X E) = — e • E. 


2 COS0B 


(23) 


(26) 
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We assume that the medium is weakly anisotropic so that 
deflection of the photon trajectory can be neglected. This 
requires that the deviation of e or /i from the unit tensor 
is small. Consider a wave propagating in the z-direction. 
The external magnetic field is assumed to be constant, 
and the medium density varies in space with character¬ 
istic length scale Hp c/u. Let E = where 

ko = Lojc. Assuming that \dA/dz\ <C ko\A\ (geometric 
optics approximation), we have 

[(eii — o,)Ax + ei2Aj, -|- eisA^], (27) 

dAy _ _ ifco _ r 

dz 2{a + 9 b) ^ 

+ (^22 — a — m sin^ 6b) Ay € 23 ^^], (28) 

A, = +62371, ^ ^29) 

€33 

where is the component of the dielectric tensor in the 
xyz coordinate system (with the z-axis along k). We shall 
hnd that near the vacuum resonance, A can vary on the 
length scale l/(fco|n+ —n_|) (see §5.2), much shorter than 
Hp, but the condition \dA/dz\ <C A:o|A| is still easily sat¬ 
isfied as long as \n+ — n_| <C 1. Using the expressions 
for tij, we find that leia/esaj ~ |(i’e — q) sin 6 *B cos 0 b | and 

_ 1 /o 

|e23/e33| — keMe SinO B / {I - Ui)\ ioT \u^ - l\ 'Z me / nip 
(and Ue ^ 1). Thus, as long as Ue <C 1 (valid for suffi¬ 
ciently low densities) and g <C 1 (valid for Bu <C 600), 
the wave is transverse for most photon energies (an excep¬ 
tion occurs when E is very close to Esi). We shall adopt 
this transverse approximation in the following Since 
a — 1, m <C 1, the evolution equation for A simplihes to 


d f Ax 
Ay 

where 


iko ail o’i2 
2 0-21 0-22 


A, \ 

Ay J 


(30) 


where An is the difference of the indices of refraction for 
the two modes and 9^ is the “mixing angle”, as given by 

—iai2 „ —2iai2 , , 

An = . , tan20^ =-. (36) 

sm20m CTii - (T 22 

When dissipation is neglected. An and 0^ are real, and 
for We 1 , we find 


An = — 


Ve cos 9b 


(tan 20 „ 


uj'^(\ — Ui) sin 20 ; 
u\A[l — Ui) sin^ 9b 


2 cos 9b 


1 - 


m 


= 13. 


(37) 


(38) 


The normal mode eigenvectors are (Ax, Ay) = 
(icos0m,sin0m) and (—isin0m,cos0m)- We dehne the 
mode amplitude A+ and A_ via 


( Ax 
\ Ay 


f i cos 0m 
V sin 0m 


+ A- 


—i sin 0 m 
cos 0m 


• (39) 


The equations for the evolution of A+ and A_ are then 
given by 



—fcoAn /2 id9m/dz f A+ \ 
—i d9m,/dz koAn/2 I A- J ’ 


(40) 

which agree with eq. (15) of LH02. Equations (35) and 
(40) have the standard forms of the level crossing problem 
in quantum mechanics (Landau 1932; Zener 1932). When 


fco |An|/2 > |(i0m/(iz| , (41) 


the polarization vector will evolve adiabatically (i.e., 
evolve along the continuous A+ or A_ curve). Evalu¬ 
ating this condition at vacuum resonance yields E ^ Ead 
[see eq. (2) and LL[02 for details]. In general, for a linear 
density profile, p = py(l — z/Hp) cz pv ex.p(—z/Hp) (for 
jzj <C Hp), the nonadiabatic jump probability is given by 


CTii = £11 — a = e' cos^ 9b + p' sin^ 0 b — a, (31) 
CT 12 = - 0-21 = £12 = i5COS0B, (32) 

022 = £22 — a — msin^ 9b = s' — a — msin^ 6b. (33) 

Equation (30) can be used to study the evolution of EM 
wave amplitude across the vacuum resonance under the 
most general conditions. Equivalently, we can use the pho¬ 
ton density matrix to study the wave propagation (§5). 


4.1. Evolution of Mode Amplitude 


It is instructive to show that when dissipation is ne¬ 
glected, eqs. (30) can be reduced to the evolution equa¬ 
tions for the mode amplitudes given in LH02. Let 


( A, 

V Ay 


= exp 


ikp 

4 


(ctii + o-22)dz 


Ax 

Ay 


(34) 


Equations (30) reduce to 


— ( 

' Ax ^ 

\ ikpAn 

cos 20m 

i sin 20m 

dz \ 

V Ay y 

1 ~ 2 

—z sin 20m 

— cos 20m 


Ax \ 
Ay J 
(35) 


.I^ump — exp 


-{tt/ 2) (E/E^df 


(42) 


Thus, in practice, when E ^ 1.3 Ead, the evolution will 
be highly adiabatic (with Ljump ^ 0.03). Figure 6 shows 
the adiabatic regime E > E^d in the B-E plane for typi¬ 
cal values of 0B and Hp. [Another way to understand the 
adiabatic condition is discussed in §5.1; see eq. (60)]. 


5. EVOLUTION OF STOKES PARAMETERS ACROSS 
VACUUM RESONANCE 

5.1. Evolution Equations for Photon Intensity Matrix 

We now study the transfer equation for the photon in¬ 
tensity matrix®, taking into account the birefringence of 
the plasma and the attenuation due to absorption and 
scattering; we relegate the derivation of the source func¬ 
tions to §6. 

The intensity matrix of the radiation is dehned as 

= C(A„A^), (43) 

where (• • •) denotes time average (a, (3 = 1,2 or x,y). 
We choose the proportional constant C such that lap 


® This approximation is also consistent with the requirement that the medium be weakly anisotropic. 

® It is beyond the scope of this paper to properly review previous works on transfer equations of polarized radiation in a magnetized medium; 
see, e.g., Unno 1956; Dolginov et al. 1970; Dolginov & Pavlov 1974; Zheleznyakov et al. 1974. These papers study the transfer equation under 
a variety of different approximations, in different contexts and with different level of generalities. 
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has the units of specific intensity. The Stokes parameters 
I,Q,U,V are related to laf^ via 


Uaff] = 


1 

2 


I+ Q U + iV 
U-iV I-Q 


(44) 


Using eq. (30), we can derive the evolution equation for 
the intensity matrix^: 

=-1^0 (^«7^7/3 + Ic.-rT+^) > (45) 

7 

where 

Taf} = -icTaf}. (46) 


Equivalently, the Stokes parameters S 
{So, Si,S 2 , S'a)''" = {I, Q, U, V)'^ evolve according to 


/ 

Q 
u 

V 

where 


dz 



- To Ti T2 To - 


( ^ \ 

hr. 

Ti To —123 ^2 


Q 

Kq 

T2 123 To —I2i 




To —122 f 2 i To 


\ V J 


(47) 


(47) 

To = —- (cTiii + cr22i) , (48) 

Tl = —- {(Tin — <T22i) , (49) 

T2 = —— {lTl2i + <T21i) = 0, (50) 

^3=2 ('^12r- ~ '721r) = <7i2r, (51) 

Vti = — {aiir — (T22r) , (52) 

^^2 = ^ (0-12r + Cr21r) = 0, (53) 

^3=2 ('^12i ~ '721i) = Cri2i. (54) 


Here aa/ 3 r = ^^{(Taf}) and <7af}i = Im((jQ,/ 3 ). Following 
Kubo & Nagata (1983), equation (47) can also be cast 
into a vector form. Define vectors in the Poincare sphere 


s = {SilSo,S 2 lSo,So/So), (55) 

T=(Ti,T2,T3), (56) 

fl = (0i,D2,03). (57) 

Then 

= /co X s + s X (T X s)]. (58) 


In the absence of damping, T = 0, and from eq. (36) we 
have® 

17 = An (cos 20^,0, sin 20m) • (59) 

With this, the adiabatic mode evolution discussed in §4 
can also be understood as follows. Consider a photon ini¬ 
tially in the plus-mode, with {Ex,Ey) oc (i cos0m> sin0m); 
it can be represented on the Poincare sphere by the vector 
s = (cos 20m, 0, sin 20m). Thus s is initially parallel to 12. 
As the photon propagates in an inhomogeneous medium, 
12 = 12/An evolves on the Poincare sphere. If 12 changes 


slowly, s will evolve in the way that it remains almost par¬ 
allel to 12 ; this requires 


d(12/An) 

dz 


< A:o|12|, 


(60) 


or 2\d9m/dz\ <C ko\^n\, in agreement with eq. (41). 

For use in §6, we define 

X^{hi,l 22 ,U,V)+. (61) 

The evolution equation can then be written as 
d 


dz 


T = -M ■ T, 


where 

[A4] = 


Clli 0 Cri 2 i /2 —CTi 2 r /2 

0 Cr22i —(Tl 2 i /2 —CTi2r/2 

— (Tl2i (Tl2i {(Tin + (722i) / 2 {(Tllr — (T22r) / “2 

— (Tl2r —(Tl2r {<T22r — (Tllr)/2 {(TinE(T22i){2 \ 


(62) 


(63) 


5.2. Numerical Results 


Figure 7 depicts four examples of the evolution of the 
Stokes parameters when monochromatic radiation (of en¬ 
ergy E = 0.1, 0.5, 1, 3 keV, respectively) propagates in 
an inhomogeneous hydrogen plasma with density profile 
P = Pv {E) e:Kp{—z/Hp), where Tip = 5 cm. The vacuum 
resonance is located at p = pv{E) and z = 0, and the 
external magnetic field is B 14 = 1 and 0 _b = 45°. In all 
cases, the radiation is assumed to be in the plus-mode and 
has intensity 1 = 1 (arbitrary units) at 2 : = —0.3 cm. For 
E = 3 keV, the evolution is highly adiabatic (see Fig. 6), 
thus the radiation will remain in the plus-mode, with Q 
changing from 1 (characteristic of 0-mode) to —I (char¬ 
acteristic of X-mode). For E = 0.1 keV, the evolution is 
highly non-adiabatic, and the four Stokes parameters re¬ 
main almost constant across the vacuum resonance. For 
E = 0.5 keV and 1 keV, we have the intermediate situa¬ 
tion where partial mode coupling/conversion takes place 
at the resonance. 

The oscillation of Q, U, V after the radiation crosses 
the vacuum resonance can be understood as follows. Since 
the modes are orthorgonal and form a complete set away 
from the resonance, we can write the radiation field after 
the resonance as 


E = 




l cos 0. 
sin 0„ 


+ A_ 


—i sin 0m 
cos 0m 


where 0m is real (i.e., the damping terms are neglected in 
calculating the modes). The corresponding Stokes param¬ 
eters are 


/=|A+p + |A_^ (65) 

Q = cos 20m (|A+p - |A_p) 

—2 sin 20m IA+II A_ I cos A/), (66) 

[/=-2|A+||A_|sinA(/), (67) 

U = sin20m(|A+|2- |H_|2) 

+2cos26*m|^+||^-| cos A^, (68) 


^ This form of the transfer equation was derived in Dolginov et al. (1970) for a general magnetoactive plasma, although no simple expression 
for Tq,^ was given there. 

^ Heyl and Shaviv (2000) used an equation similar to eq. (58) (with T = 0) to study the polarization evolution of EM waves in a pure 
magnetized vacuum (for which dm = 0) with a nonuniform external magnetic field. 
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where |^±| oc exp(—/9 k±z/ 2) is the mode amplitude (k± 
is the opacity) and A(j) = A^o + ko J (n+ — n-)dz is 
the phase difference between the modes (A(/)o is the ini¬ 
tial phase shift). Clearly, Q, U, V oscillate on the length 
scale ~ l/[fco|^-i- — n-\]. Also note that away from the 
resonance, 0m ^ 0 or tt/2, so that the oscillation of Q 
diminishes rapidly after the resonance. 

The results of Fig. 7 are obtained using eq. (62). We 
can obtain the same results by integrating eq. (30) and 
calculating the Stokes parameters from eq. (43); this re¬ 
quires many more integration steps than using eq. (62) 
directly to achieve the same numerical accuracy. Since 
the photon modes are well-defined away from the vacuum 
resonance, we can calculate the mode amplitudes A_|_ and 
A- from Ax and Ay [from integration of eq. (30)] us¬ 
ing eq. (64). Figure 8 shows the evolution of the ratio 
|A+|/|A|, where |A| = (|A+p -I- |A_p)^/^. Using eq. (42), 
we find that Pjump = 1 — 3.2 x 10“^, 0.84, 0.27, 6.5 x 10“^ 
for E = 0.1, 0.5, 1, 3 keV, respectively. Starting with 
|A_|_|/|A| = 1 before resonance, this implies |A+|/|A| = 
(1 — Fjump)^^^ = 5.7 X 10“^, 0.40, 0.85, 1 after the reso¬ 
nance, in agreement with the numerical results shown in 
Fig. 8. Note that except for the E = 0.1 keV case, the 
three other cases all he in the large Faraday depolariza¬ 
tion regime (see Figs. 4, 5), and the agreement between 
Fig. 8 and eq. (42) is not surprising. 

Figure 9 shows another example of the evolution of po¬ 
larized radiation, with E = 1.7 keV, Bu = 5, 0b = 45° 
and E[p = 5 cm. For this case, [3 = 1.066 1 at p = pv, thus 
the two modes nearly collapse at the vacuum resonance. 
Nevertheless, we find that eq. (42) provides an accurate 
description for the evolution across the resonance: after 
the resonance, |A_| = Pj^mp = 0-93 and decays slowly, 
while the plus-mode damps quickly because of the large 
opacity. 

We note the distance over which resonant evolution oc¬ 
curs is much smaller than the the photon mean free path 
(see Figs. 7-9). We can estimate the width of the resonance 
as follows. From eq. (38) we have 

5- n - 1 - PV 9) ■ 69) 

2 cos0B 

At the resonance (3 = 0. If we dehne the edge of the reso¬ 
nance region using the consition \ j3\ = (3ed k, a few, we find 
the half-width of the resonance to be given by (LH02) 
2(3edHp cos 0B 

= -TJ27, - 

Ue 11 - Ui I sm 0B 


:xl.73x 10-^(3ed\l-u^\-^^^^ f^)/7p.(70) 

sim 0B \Bi4/ 

On the other hand, the imaginary parts of the indices of 
refraction of the two modes are given by 

Im(n+) ~ 7 eUe sin^ 0 b cos^ 0^, (71) 

Im(n?.) ~ 7 eUe sin^ 0 b sin^ 9m, (72) 


[valid for Ue ^ 3 and Uj <C 1; cf. eq. (25)1. The mean 
free path I is given by l~^ = (a;/c)Im(n^), and thus 
I > c({ujjeVe) (the right-hand side being the mean free 
path at B = 0). When collisional damping dominates, 
7 e ~ 7ei » Jre [see eqs. (lO)-(ll)], we have 


,> 0 26__ 

-zvi(i - 


cm. 


(73) 


When scattering dominates, I > 2.5pi ^ cm. Thus the re¬ 
lation I » {Az)v is almost always satisfied. 

6. TRANSFER EQUATIONS FOR STOKES PARAMETERS 

Including the source term, the transfer equation reads 
[cf. eq. (62)] 

^X = k.VX=-M-X + 5em + 5,c, (74) 

as 

where we have used ds (instead of dz) to specify the deriva¬ 
tive along the ray. We derive the source functions and 
Ssc below. 

6.1. Source Function: Emission 


The emissivity <Sem can be derived by considering the 
radiation field in thermodynamic equilibrium with matter, 
for which X = (i?,>/2, i3i,/2, 0,0)+. Note that under this 
situation, Ssc should completely cancel out the scattering 
contribution to M. ■ X (see §6.2.2). Thus 

Sem = (cnu, 0 - 222 ,0, -2CTl2r)^,„. (75) 

This is equivalent to adding to the RHS of eq. (47) a source 
term 

-koB4To,Ti,T2,n)-^^. (76) 

This is also the same as adding to the RHS of eq. (45) a 
term 

1 


■^koB^ 


f 'T' ^ 

yj. cxp “T X J 


(77) 


The subscript “em” in eqs. (75)-(77) implies that the terms 
proportional to jre = Iri should be excluded in evaluating 
these equations, since these terms are related to scattering 
contributions (§6.2.2). Equations (31)-(33) yield 
0-1 li = £i cos^ 0B + Pi sin^ 9b, 


0-222 — Si, 

0^12r — 9i^S‘^9B. 

From eqs. (7) and (8), we Hnd 
£i±gi = Im(£ ± 5 ) = 


A± 


Ve'yfi + (1 T Uy^)^Ue7re + (1 ± ul^'^)'^Vilr 


-1 


(78) 

(79) 

(80) 


,(81) 

(82) 

(83) 

(84) 


— 


kpVeBy 


Pi = Im(r7) = Ue7ei + Ue7e + 
where 

A±=[(l±„y 7 (lT»‘'U+ 7 t 

and 

7± = it + (1 ± + (1 T 

Thus eq. (75) becomes 

cos20b (A+7+ -f A_7“) -h 27ei sin^0B 

^+lti + A-7e“2 

0 

2 (A+ 7 J) - A-7e“2) COS0B 

(85) 

Note that the terms in eq. (85) can be easily related 
to the free-free opacity: in magnetic fields, a photon of a 
certain polarization has an absorption opacity given by 
K® = K+Je+l^-b K-le_l^-t Kojeo]^, (86) 

where eo = ez, e± = {ex ± iev)!'/^ are the spherical 
components (with B along the Z-axis) of the photon’s unit 
polarization vector e, and 

^0 ± . kp , , 

K± = —VelmA±, Ko = —Velei- (87) 

Equation (87) agrees with (and generalizes to include the 
Z > 1 case) the result of Potekhin & Chabrier (2002). 



6.2. Source Function: Scattering 

We first derive the source function due to electron scat¬ 
tering. The contribution from ion scattering can be easily 
added (see below). 

Consider an incident EM wave with electric field 

E'(r,t) = where E' = + E'yey,. In 

the coordinate system XYZ with B along the Z-axis, the 
direction of the wave vector k' = fepk' is specified by the 
polar angles 9'^ and (j)'g. We define the directions of the 
unit vectors e^' and e^/ such that lies in the merid¬ 
ian plane and points in the direction of increasing 9'^, and 
e^' X Byf = k'. Similarly, the scattered wave is denoted 

by E(r, t) = E where E = + EySy, and the 

polar angles for k are 9b and (j)B- We also define unit 
vectors 

e± = (ejf ± ier)/V2, eg = ez- (88) 


We can write the incident wave as 


E' = 


E'e*, with El =e„-E'. 


(89) 


a=0,± 

(Here and henceforth, refers to sum over a = 0, ±). 
The scattered wave is given by the standard dipole formula 
E = k X (k X d)/(c^r), with 


d = c^re'V(l-I-OM^^) ^EF 


(90) 


where re = ^ j(me<?\ Thus the Jones matrix relating E' 
and E is given by 


Ex 

E -,1 


h 

r 


a b 
c d 


EL 

E'y 


where 


a = ^(l + Q!uy^) ^{ex-e*^){ea-ex'), 

Oc 

b = '^{l + aul/'^)-^{ex-e*J{ea • e^/), 

Oc 

c = ^(1 -f aul^^y^ey ■ e;)(e„ • e„/), 

Oc 

d = ^(1 -h aul/^)-\ey ■ e*)(ee, ■ ey^). 


(91) 

(92) 

(93) 

(94) 

(95) 


Using Bx = cos 9 b( cos (j^B^x + sin^BOy) — sin^sez and 
By = — sin(/>Bex + cos</>bGy (and similarly for Bx' and 
by/), we find 

a = cos 9'g cos 9 b + sin9'gSm9B, (96) 


b = cos 9B ( A(j }), 

c = — cos 9'b (A(j )), 

d = C^^\A<p), 


where® 


C(®)(A</)) = i 


eiA<j> 


^ — iAcp 


2 \l-ul/^ 1 


S<->(A« = 1 ( 


oiA(f} 


, 1/2 
f Me 


(97) 

(98) 

(99) 

( 100 ) 




with = (pB — 4 >'b- Using the definition of X [see 
eq. (61)], we can relate X(k) for the scattered radiation 
to I' — X(k') for the incident radiation. Integrating over 
all possible incident directions k', we obtain the source 
function due to electron scattering: 

Sif (k) = uerl jdVL' R(®) (k' ^ k) • X(k'), (102) 

where dVL' = d^k' = dcos9'gd(j)'g, and the scattering phase 
matrix R*-®^ is given by 
R(®)(k' ^ k) = 

|ap l&P (a*b)r (a*b)i 

|c|® jdp (c*d)r (c*d)i 

2 (a*c)r 2(b*d)r (a*d + b*c)r (a*d—b*c)i 

_—2(a*c)i —2(b*d)i —(a*d + b*c)i (a*d — b*c)r 

In eq. (103), we have used the notation (a*b)r = Re(a*6) 
and (a*b)i = Im(a*&), etc. 

The source function due to photon-ion scattering can be 
similarly written as 

5«(k) = Jdn' R(*)(k' ^ k) • X(k'), (104) 

where is the ion number density, = (Ze)"^/(AnipC^), 
and R*^*^ has the same form as eq. (103), except that in the 
expressions for a,b,c,d [eqs. (96)-(99)] one should replace 
C^^\A(j)) and S^^\A(j)) by 


(103) 


C«(A</)) = - 


S^^^AP) = 






1 -I- M, 


1/2 


1 — M, 


1/2 


oiAcp 


V 1 + 1 - u) 

It is easy to check the source functions derived above re¬ 
duce to the appropriate zero-field results (Chandrasekhar 
1960) when R —> 0. 

6.2.1. Special Case: Axisymmetry 

When X is independent ol pB [as is the case when B 
is perpendicular to the (local) stellar surface], the source 
function reduces to 

5(®/k) = 27rner® Jdpc' R(®)(k' ^ k) • X(k'), (107) 

where fj! = cos9'g. The azimuthal average of R*^®) is given 
by 

(Ja]®) yc^Fe 0 yL'ye^G 

fi'^Fe Fe 0 fl'Ge 

0 0 0 0 

2 p'^y,Ge 2/iGe 0 2 p'yLFe, 

where = cos 9b and 

{\a\^) = (f,'f,)^Fe + (l-fi'^)(l-p^). 

The functions Fg and Ge in eqs. (108)-(109) are 

Fe = (1G(®)(A</))1®) = (1R(®/A0)1®) = ^ + ^ 


^—iA(f) 

JG 


(105) 


(106) 


R(®/k' ^ k) = 


® To avoid divergence at Ue = 1, one can replace (1 + caAJ'^') by (1 + 

1 /2 

relevant equations. Similarly, one can replace (1 + ocu^ ) by (1 + 27 ^ 
these damping terms explicitly in the equations of this subsection. 


^ , ( 101 ) Ge = f( g(®) (AP) [r(®)(Ac/.)] 


2(1-Me)2 

„l/2 


(108) 

(109) 

a 10) 


n 12- (111) 

1/2 (l-«e)2 

27 e + oluJ ) (where 7 e is the dimensionless electron damping rate) in all 
1 /2 

+ olu^ ) to avoid divergence at -Uj = 1. For simplicity, we do not include 
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The source function due to photon-ion scattering is 

5«(k) = 27Tn,rf Jdfi’ ^ k) • X(k'), (112) 


where has the same form as eq. (108), except that the 
functions Fg and Gg are replaced by 


1 + Ui 

2 {l-UiY' 


,1/2 


G. = 


(l-w,)2- 


(113) 


6.2.2. Special Case: Thermodynamic Equilibrium 


When the radiation field is in complete thermodynamic 
equilibrium with matter, X = {B^I2, 5,^/2, 0 , 0 )+, and the 
source functions due to electron and ion scatterings be¬ 
come 

r^^Fg + {l- p^)/2- 


Si'g>{k)=ngaTB, 


5«(k) 


Fg 

0 

2 piGg 

fx^F, + (1 - p^)/2 

F^ 

0 

2 plG, 


(114) 


(115) 


where ar = 87 rrg /3 and cr^ = Sirr^/S. On the other 
hand, the sink term due to scattering is 

^0 


iM-Xl^ = 


(cili, (7221, 0,-2ai2r)l, (116) 


where the subscript “sc” implies that one should keep 
only the terms proportional to jre and jri in eqs. (81)- 
(82) (the terms proportional to jgi are due to absorp¬ 
tion). Neglecting the damping 7 ^ in eq. (83), we find 
e* = 2{vgjreFg + Vl'yriFi), Pi = 2{vg-frgGg + VijriGi), and 
Vi = Vejre + Vilri- Substituting these into eqs. (78)-(80) 
and using 


Vgjre = ngar/ko, Vijri = riialp'’/ko, (117) 

we can easily show that — (A4 • = 0 is 

satished, as it should. 


7. CONCLUSION 

The main results of this paper can be summarized as 
follows. 

1. We presented a physical discussion of the effect of vac¬ 
uum polarization on radiative transfer in strong magnetic 
fields and thermal radiation from magnetars (§2). This 
substantiated and clarified our previous results (LH02 and 
HL03): vacuum polarization depletes the high-energy tail 
of the thermal spectrum and suppresses spectral lines. 

2. We studied the mode properties and examined sev¬ 
eral issues associated with the resonance due to vacuum 
polarization in a magnetized electron-ion plasma (§3). We 
showed that for superstrong magnetic fields, mode collapse 
and the breakdown of Faraday depolarization occur near 
the vacuum resonance for a wide range of propagation di¬ 
rections. Therefore a rigorous treatment of radiative trans¬ 
fer in magnetar atmospheres must go beyond the modal 
description of the radiation. 

3. We studied the propagation of polarized radiation in 
the inhomogeneous atmosphere of magnetars (§4 and §5). 
We found that even with mode collapse and breakdown of 
Faraday depolarization, the evolution of polarized X-rays 
across the vacuum resonance can be described by a jump 
probability from one mode into another. 

4. We derived the general transfer equations for the 
photon intensity matrix (Stokes parameters) in the mag¬ 
netized plasma-vacuum medium (§ 6 ). Explicit expressions 
for the source functions due to emission and scatterings 
were obtained. These transfer equations are useful for ac¬ 
curate calculations of radiation/polarization spectra from 
magnetars. 


We thank G. Pavlov and A. Potekhin for useful dis¬ 
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remarks. This work is supported in part by NASA grants 
NAG 5-8484 and NAG 5-12034, NSF grant AST 9986740, 
and a fellowship from the A.P. Sloan foundation. 
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Fig. 1. — A schematic diagram illustrating mode conversion due to vacuum polarization: As a minus-mode (plus-mode) photon, which 
manifests as X-mode (O-mode) at high density (for E > Ebi), traverses the vacuum resonance density py, it will stay as minus-mode 
(plus-mode) and become O-mode (X-mode) at low density if the adiabatic condition is satisfied (see also the left panels of Fig. 3). In this 
evolution, the polarization ellipse rotates 90° and the photon opacity changes significantly. 
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Fig. 2.— A schematic diagram illustrating how vacuum polarization-induced mode conversion affects radiative transfer in a magnetar 
atmosphere. When the vacuum polarization effect is turned off, the X-mode photosphere (where optical depth 1) lies deeper than the 
O-mode. With the vacuum polarization effect included, the X-mode effectively decouples (emerges) from the atmosphere at the vacuum 
resonance, which lies at a lower density than the (original) X-mode photosphere. 
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Fig. 3. — The polarization ellipticities K [upper panels; see eq. (17)] and refractive indices (lower panels) of the photon modes as functions 
of density near vacuum resonance for S = 5 x 10^"^ G, 6b = 45°, and Ye = 1. The left panels are for E = 5 keV, and the right panels 
for E = 1 keV. On the right panels, the light curves show the results when damping is neglected, while the heavy lines include damping 
(we set T = 5 X 10® K in evaluating 7ei)- On the left panels, the results are indistinguishable with and without damping. Note that the 
correspondence between the plus-, minus-modes and the X-, O-modes depends on the density and the photon energy (E > Ebx vs. E < Ebi)- 
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Fig. 4.— The critical angle 6b = ^coll for mode collapse at vacuum resonance as a function of photon energy. The four curves correspond 
to Bi 4 = 0.1, 1, 5, 10. Note that for E = Ebi (the ion cyclotron energy), ^coii is very close to 90^^. For ^coii < 6b < 180° — ^coii? the modes 
are highly nonorthogonal and the condition for Faraday depolarization breaks down (see also Fig. 5). 
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Fig. 5.— The critical photon energy E = Ob) for the breakdown of Faraday depolarization. For E ^ F’rar (the unshaded region), 

the condition for large Faraday depolarization is satisfied. Note that the boundary E = F'far(-B, Ob) also corresponds to Ob = ^coll(-^> (see 
Fig. 4). 
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Fig. 6.— The critical photon energy for adiabatic mode conversion as a function of magnetic field strength. The unshaded region refers to 
the adiabatic regime in which the condition E > is satisfied [see eq. (2)]. The solid line is for Oq = 45° and the dashed line for Oq = 15°, 
both assuming Hp = 5 cm. 
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Fig. 7.— Evolution of Stokes parameters across the vacuum resonance. The parameters are B 14 = 1 , = 45°, and the density profile is 

p = py ex.p{—z/Hp) with Hp = 5 cm. At 2 = —0.3 cm, the radiation is in the plus-mode, with / = 1. The solid lines are for E = 1 keV, the 
short-dashed lines for E = S keV, the dotted lines for E = 0.5 keV, and the long-dashed lines for E = 0.1 keV. 
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Fig. 8 .— Evolution of the ratio |A+|/|A| for the cases depicted in Fig. 7. Here \ A±\ is the amplitude of the plus and minus-modes and 
A| = V|A+P + !A_p. 
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Fig. 9.— Evolution of Stokes parameters J, Q and mode amplitudes |A+|, |A_| across vacuum resonance, for E = 1.7 keV, B 14 = 5, 
Oq = 45°, and Hp = 5 cm. At 2 = —0.03 cm, the radiation is in the plus-mode, with A^. = 1. The solid lines show the results including 
realistic damping terms in the dielectric tensor, and the dotted lines show the results when the damping terms are turned off. Note that for 
this case, (3 = 1.066 2 at p = py ^ so that the mode amplitudes are well defined only away from the the resonance (A+ and A- are calculated 
assuming that 9m is real and the two modes are completely orthogonal). 





















